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Abstract. The equation for QCD string proposed earlier is reviewed. This equation 
appears when we examine the gonihedric string model and the corresponding transfer 
matrix. Arguing that string equation should have a generalized Dirac form we found 
the corresponding infinite-dimensional gamma matrices as a symmetric solution of the 
Majorana commutation relations. The generalized gamma matrices are anticommuting 
and guarantee unitarity of the theory at all orders of v/c. In the second quantized form 
the equation does not have unwanted ghost states in Fock space. In the absence of 
Casimir mass terms the spectrum reminds hydrogen exitations. On every mass level 
r = 2,4, .. there are different charged particles with spin running from j — 1/2 up to 
jmax — r — 1/2, and the degeneracy is equal to d r = 2? — 1 = 2j max . This is in contrast 
with the exponential degeneracy in superstring theory. 

Introduction. There is some experimental and theoretical evidence for the 
existence of a string theory in four dimensions which may describe strong inter- 
actions and represent the solution of QCD [1]. 

One of the possible candidates for that purpose is the gonihedric string which 
has been defined as a model of random surfaces with an action which is propor- 
tional to the linear size of the surface [4] 

A{M) = m J2 A u ' °( a ) = k - °l c . (!) 

<ij> 

where Ay is the length of the edge < ij > of the triangulated surface M and 
Oij is the dihedral angle between two neighbouring triangles of M sharing a 
common edge < ij >. The model has a number of properties which make it very 
close to the Feynman path integral for a point-like relativistic particle. In the 
limit in which the surface degenerates into a single world line the action becomes 
proportional to the length of the path and the classical equation of motion for 
the gonihedric string is reduced to the classical equation of motion for a free 
relativistic particle. At the classical level the string tension is equal to zero and, 
as it was demonstrated in [4] , quantum fluctuations generate the nonzero string 
tension o quantum = ^ (1 — ln4), where d is the dimension of the spacetime, /3 
is the coupling constant, a is the scaling parameter and <; = (d — 2)/d in (1). 

It is natural therefore to ask what type of equation may describe this string 
theory in the continuum limit. The aim of the articles [5, 6] was to suggest 
a possible answer to this question. The analysis of the transfer matrix shows 
[5] that the desired equation should describe propagation of fermionic degrees of 
freedom distributed over the space contour. When this contour shrinks to a point, 
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the equation should describe propagation of a free Dirac fermion. Thus each 
particle in this theory should be viewed as a state of a complex fermionic system 
and the system should have a point-particle limit when there is no excitation of 
the internal motion. In the given case this restriction should be understood as 
a principle according to which the infinite sequence of particles should contain 
the spin one-half fermion and the equation should have the Dirac form [5, 6] 



{ i r„ &* - M } & =0. (2) 

The invariance of this equation under Lorentz transformations x'^ — x v , 
0{A) \P(x) leads to the equation for the gamma matrices [8, 9] F v = AJ 1 0^ 1 . 
If we use the infinitesimal form of Lorentz transformations A^ v = rj^+e^, = 
1 + 2 € t^ ^ v ^ follows that gamma matrices should satisfy the Majorana com- 
mutation relation [9] 

[r^, hp] = ?? m a r p - rjpp A, (3) 

where I pv are the generators of the Lorentz algebra. These equations allow us 
to find the matrices when the representation of the 1^ is given. 

Ettore Majorana in 1932 [9] found a solution of the equations (3) which 
differs from the Dirac gamma matrices [8]. The original Majorana solution for 

matrices is infinite-dimensional (see equation (14) in [9]) and is given by the 
formulas: 

< j, m\ r Q \j, m >= j + 1/2 
<j,m\ r z |j + l,m>= ^y/(j + m + l)(j-m + l) 

< 3,m\ r z \j - l,m >= - % -yJ (j + m)(j - m) 



<j,m\ r+ |i + l,m-l>= -y/{j-m + l){j-m + 2) 
< j,m\ r + \j - l,m- 1 >= l -^J{j + m)(j + m - 1) 



< j,m\ T_ \j + l,m+ 1 >= -^VU + m + + m + 2) 



< j,m\ -T- \j - l,m- 1 >= -m)(j - m-1). 

One can see from this solution that the mass spectrum of the theory is equal to 

- wm- <4 » 

where j = 1/2,3/2,5/2, .... in the fermion case and j = 0, 1,2, .... in the boson 
case. The main problems of Majorana theory are the decreasing mass spectrum 
(4) , absence of antiparticles and troublesome tachyonic solutions - the problems 
common to high spin theories [10, 16]. 
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An alternative way to incorporate the internal motion into the Dirac equation 
was suggested by Pierre Ramond in 1971 [11]. In his extension of the Dirac equa- 
tion the internal motion is incorporated through the construction of operator- 
valued gamma matrices. The equations which define the matrices are 

< W >= 7 „, {r„(r), r„(r')} - 2v, S( ^-j (r - r) ), 

r+(r) 7o = 7o r^r), where it is required that the proper-time average < ... > 
over the periodic internal motion with period 27ra = 1/a should coincide with 
the Dirac matrices. The solution has the form 

oo 

r n( T ) = 7m + H5 ^2 [b™exp(-inu)T) + b+ n exp(inoT% 

71=1 

where 6's are operators obeying the anticommutation relations {b™,b+ m } = 
^g^uSnm and all others are set equal to zero (co = 1/a). The mass spectrum lies 
on linear trajectories and the free Ramond string is a consistent theory in ten 
dimensions and the spectrum contains a massless ground state [2, 3]. 

In both cases one can see an effective extension of Dirac gamma matrices into 
the infinite-dimensional case, but these extensions are quiet different [5, 6]. For 
our purposes we shall follow Majorana's approach to incorporate the internal 
motion in the form of an infinite-dimensional wave equation. Indeed in [5, 6] 
the Majorana theory was interpreted as a natural way to incorporate additional 
degrees of freedom into the relativistic Dirac equation. Unlike Majorana the 
authors consider the infinite sequence of high-dimensional representations of the 
Lorentz group with nonzero Casimir operators (a • b) and (a 2 — b 2 ). These 
representations (jo, A) and their adjoint (jo,— A) are enumerated by the index 
r = jo + 1/2, where r = 1, ...,N and jo = 1/2,3/2, ... is the lower spin in the 
representation (j , A), thus j = jo, jo + 1, We took the free complex parameter 
A in the real interval —3/2 < A < 3/2 in order to have real matrix elements 
for the Lorentz boost operator b. These representations are infinite-dimensional 
except of the case jo = 1/2, A = ±3/2. 

The duality transformation (jo! A) — > (A; jo), defined in [5], leads to a 
subsequent restriction on the free parameter A and requires A = 1/2, so that 
the dual representations <9 r , Of become finite-dimensional (1/2, ±(1/2 + r)). The 
corresponding equation found in [5] is not in contradiction with no-go theorem 
of [16], because dual representations are finite-dimensional. However having a 
physically acceptable spectrum this equation admits unitarity only at zero order 
ofv/c [5]. In the article [6] we found a new equation which corresponds to a sym- 
metric solution of the Majorana commutation relations and admits unitarity in 
all orders of v/c. It is based on the same dual representations (1/2, ±(1/2 + r)) 
of the Lorentz algebra and is a natural extension of the previous equation [5] 
because the new gamma matrices have the same form as the old ones but with 
additional antidiagonal elements. The new gamma matrices are anticommuting 
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and satisfy generalized Clifford algebra [6] 

{r^r u } = 2 Vtlv r 2 . (5) 

They guarantee unitarity of the theory at all orders of v /c. The anticommutation 
relations (5) are covariant, because the matrices F = A" F v satisfy the same 
relations (5) as one can check by direct computation (notice that d-T^ — r 2 ). 
For the completeness we shall review the logical and analytical steps which 
lead to equation presented in [5] and then derive the new equation which has 
anticommuting gamma matrices [6]. 

Lorentz algebra. In terms of 5*0(3) generators a and Lorentz boosts b (a x = 
il 23 a v = il 31 a z = il 12 b x = il w b v = il 20 b z = il 30 ) the algebra of the 
5*0(3, 1) generators can be rewritten as [9, 15, 8] (we use Majorana's notations) 

[a x ,a y ] = ia z [a x , b y ] = ib z [b x , b y ] = -ia z . (6) 

The irreducible representations of the 50(3) subalgebra (6) are well known. 
The dimension of is 2j + 1 and j = 0, 1/2, 1, 3/2, ... The representation 

= (jo', A) of the Lorentz algebra can be parameterized as the sum [15, 9, 8] 
0(jo, A) = ©X^jlj where jo defines the lower spin in the representation 
and A is a free complex parameter. The A appears as an essential dynamical 
parameter which cannot be determined solely from the kinematics of the Lorentz 
group 1 . The adjoint representation is denned as — (jo; —A). We shall consider 
the case r = (r — 1/2, A) and —3/2 < A < 3/2 to have real matrix elements of 
the operator bk, k — x,y,z [5]. 

Superposition principle and physical constraints. The Majorana commu- 
tation relation (3) together with the last equations allow to find matrices 
when a representation of the Lorentz algebra is given [9]. Because Jo 
commutes with spatial rotations a (see (3)) it should have the form 

< j, m| r<r' \j'm' >= jf ■ Sjj, ■ 8 mm , r, r' = N, i, 1, N (7) 

where we consider N pairs of adjoint representations — (Ojy, ■■■ ,0^,01, - ■■ , 0n) 
with j = 1/2, ...,N ~ 1/2. Thus j rr ' is 2N x 2N matrix which should satisfy 
the equation for Jb (3) [9] 

r b 2 z - 2 b z r b z + b 2 z r Q = -r . (8) 

If and are two solutions of the equation (8), then their sum is also a 
solution of (8) [5] and we shall use this superposition property in order to find 

1 The representation is finite-dimensional if A = jo + r, r = 1, 2, 3, .... The representa- 
tions used in the Dirac equation are (1/2, —3/2) and (1/2, 3/2) and in the Majorana 
equation they are (0, 1/2) in the boson case and (1/2,0) in the fermion case. The 
infinite-dimensional Majorana representation (1/2,0) contains j = 1/2,3/2,... mul- 
tiplets of the SO(3) while (0, 1/2) contains j = 0, 1, 2, ... multiplets. 



The QCD string and generalized wave equation 



5 



solution with required physical properties. In [5] the authors were searching the 
solution of the above equation in the form of Jacoby matrices 



NN-l 



7j 



\ 



,N-1N 



„N-lN-2 



N-1N-2 



V 





,NN- 



N-1N 



1i 





*4 



J 



N 
j. 
N-l 



1> 



(9) 

In the subsequent work [6] we found the finale solution which has additional 
nonvanishing antidiagonal elements 7j r and represents a symmetric solution of 
the equation (8) of non-Jacobian form. 

These solutions of the equation (8) are defined up to a set of constant factors 
which are independent from j. Indeed, because Jacoby matrices (9) have a spe- 
cific form, the original equation (8) factorizes into separate equations for every 
element Y j r+1 of the Jacoby matrix and the solution has the form [5] 



7. 



rr+1 



Const 



3 + 3 
4r 2 - 1 



(10) 



It has (4iV — 2) j-independent free constant (this fact reflects the superposition 
property of the equation (8)). This freedom allows us to impose necessary phys- 
ical constraints on a solution requiring [5, 6]: 

— i) physical behaviour of the spectrum, 

— ii) Hermitean property of the system, 

— Hi) reality and positivity of the probabiity density matrix 

This means that we require sensible physical behaviour of the theory following 
one-particle interpretation of the solutions [8, 12]. 

In order to impose these constraints one should study the spectral properties 
of the matrices of infinite size with matrix elements 7j r+1 and 7", which have 
complicated "root" dependence. The first inspection of the solution (10) simply 
shows that every element 7j r+1 grows like w j and in general all eigenvalues 
Cj will also grow with j. Therefore the mass spectrum Mj = M/ej will have 
Majorana-like behaviour (6) Mj w M/j. To avoid this general behaviour of the 
spectrum one should carefully inspect eigenvalues of the matrix i~o for small 
values of N and then for arbitrary N [5]. The parameter N plays the role of a 
natural regularization. The B — H — £ — Si — ^-solutions which appear (see [5] 
and below) have exeptional behaviour: half of the eigenvalues of the spectrum 
are increasing and the other half are decreasing. One can achieve this exeptional 
behaviour of the solution by tuning the free constants in the general solution 
(10). However solutions B — H — S — Si — U 2 cannot be accepted [5] because 
still half of the eigenvalues produce a mass spectrum which has an accumulation 
point at zero mass. 
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Dual representations. Last phenomenon can be understood on the example 
of the Dirac equation. For that let us define the dual representation as [5] = 
(jo; A) — ► (A; j ) = dual . This symmetry transformation imposes constraints 
on the free parameter A, so that it should be integer or half-integer. This is 
because in the representation (j ; A) j must be integer or half-integer [14]. For 
the dual representations and dual the matrix elements of Lorentz generators 
Ipu are precisely the same, the only difference between them is that the lower 
spin is equal to jo for the reperesentation and is equal to A for its dual dual . 
Therefore any solution of the Majorana commutation relations (3) for can 
be translated into the corresponding solution r dual for dual by exchanging j 
for A [5]. 

The dual transformation of the Dirac representations (1/2, —3/2) and (1/2, 3/2) 
would be infinite-dimensional (3/2,-1/2) and (3/2,1/2) with j = 3/2,5/2,... 
and the corresponding solution J^ uoJ has the form jj 1 = -fj 1 = j + 1/2 with 
the following mass spectrum 

M ?ir«c *«u = _M_ j = 3/2,5/2,.... (11) 
3 + !/2 

This Majorana-like mass spectrum is dual to the physical spectrum of the Dirac 
equation 

Mf irac = M, .7 = 1/2. (12) 

The dual equation is simply unphysical, but we have to admit that the whole de- 
creasing mass spectrum of the dual equation corresponds or is dual to a physical 
Dirac fcrmion. From this point of view we have to ask about physical properties 
of the equations which are dual to "unphysical" ones B — H — £ — Si — S 2 . 
The dual transformation completely improves the decreasing mass spectrum of 
these equations [5], as it takes place in (11) and (12), and indeed the last ZJ dual - 
equation has the spectrum of particles and antiparticles of increasing half-integer 
spin lying on quasilincar trajectories of different slope. However having a physi- 
cally acceptable spectrum (see constraint i) ) this equation admits unitarity (see 
constraint ii) ) only at zero order of v/c. The equation which admits unitarity 
at all orders have been found in [6] and we shall review it below. 



Invariant product and probability density. Let us introduce the invariant 
scalar product < &i \ W 2 > = < | #2 >, where = I + \t^ v and 
the matrix [2 is defined as < ^ | <F 2 >= &i &2 = &f O &2 = J m ^ ., m , W 2 f 
with the properties 

n a k = a k n nb k = b+ n n = q + . (13) 

From the first relation it follows that fl = uj rr ■ Sjj> ■ 8 mm > and from the last two 
equations, for our choice of the representation and for a real A in the interval 
—3/2 < A < 3/2, that wj r = uj rr = 1 w| = 1, thus £1 is an antidiagonal 
matrix. The conserved current density is equal to J M = <P 7^ J p = 0. 
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The current density J should be real and positive definite (see constraint iii) ) , 
which is equivalent to the requirement that 

r+Q = nr^ (w) 

and to the positivity of the eigenvalues of the matrix p = J? Jo. 

The B — H — £ — S\ — ^-solutions. The basic solution ( B-solution) of 
the equation (8) for the Jo has the form (9), (10) with all set of constant factors 
equal to Const = i [5] and jj 1 = -fj 1 = j + 1/2. The positive eigenvalues ej 
can be found in [5] . They show that the coefficient of proportionality behind j 
drops N times compared with the one in the Majorana solution ej = j + 1/2 in 
(6) and now many eigenvalues are less than unity and the corresponding masses 
Mj = M/ej arc bigger than the ground state mass M. This actually means 
that increasing the number of representations in = (0 N , • • • , Q\, ©i, • • • , On) 
one can slow down the growth of the eigenvalues. To have the mass spectrum 
bounded from below one should have spectrum with all eigenvalues Cj less than 
unity. 

In the limit N — > oo the B-solution is being reduced to the form 

7} +1 r = ^ r+1 = -rr 1 * = ^ r+1 = 1 \f^i~b -^ + 1 /2 (15) 

and 7j 1 = 7* 1 = j + 1/2, where r = 1, 2, .... All eigenvalues Cj tend to unity 
when the number of representations N — > 00. The characteristic equation which 
is satisfied by the gamma matrix in this limit is 

id) - 1) J+1/2 = . 7 = l/2, 3/2, ,5/2,-- • (16) 

with all eigenvalues ej = ±1. Therefore all states have equal masses Mj = 1, 
but the Hamiltonian is not Hermitian (_Tq~ 7^ 7q) 2 ■ The matrix Q Jo has the 
characteristic equation (w^ 7j — l) 2 i +1 = with all eigenvalues equal to pj = +1. 
Thus the matrix ft r is positive definite and all its eigenvalues are equal to one, 
but the relations Q To ^ Jq" i7, Jq" 7^ Jo do not hold. What is crucial here is 
that we can improve the B-solution without disturbing its determinant which is 
equal to one (16) (Detr n = 1). The last property of the determinat is necessary 
to keep, in order that the spectrum will be symmertically distributed arround 
unity. 

The Hermitian solution (H-solution) of (8) for Jo can be found as a phase 
modification of the basic B-solution [5] 



-) 3 > r+1/2. (17) 
Tr 7^ = 2.7 + 1, thus e? • 



r+l r _ _ r r+1 _ _ r+1 r _ r r+1 _ ■ 1( 3 + 3 



The determinant and the trace are equal to Det 7j = ±1, 
- ' S 2 +i/2 = 1. e ? + - + S 2 +i/2 =3 + V2. 
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These matrices are Hermitian = Jb, but the characteristic equations are 
more complicated now. These polynomials p(e) have reflective symmetry and 
are evenpj(e) = e 2j+1 pj(l/e), Pj(—e) = Pj(e) therefore if ej is a solution then 
1/ej, — ej and —1/cj are also solutions 3 . The eigenvalues ej can be found in 
[5]. The changes in the phases of the matrix elements (17) result in different 
behaviour of eigenvalues. The half of the eigenvalues (decreasing eigenvalues) 
produce quasilinear trajectories with nonzero string tension and the other half 
(increasing eigenvalues) affect the low spin states on trajectories, so that the 
smallest mass on a given trajectory tends to zero (see [5]). The matrix ft Jo 
has again the characteristic equation (uij -fj — 1) 2: > +1 = and all eigenvalues 
are equal to one. Thus again the matrix ft Jo is positive definite because all 
eigenvalues are equal to one, but the important relation ft Jo 7^ r£ ft does not 
hold. 

The solution of (8) for Jp with both properties -T + = Jo and ft Jo = r o + ft 
can be found by using the basic solution rewritten with arbitrary phases of the 
matrix elements and then by requiring that Jo should be Hermitian = Jo 
and should satisfy the relation ft Jo = r o + ft. This solution, U-solution, is 
symmetric and has the form [5] 

7} +lr = ^ r+1 =^ +lr = ^ r+1 = \[(£S[-^ J>r+l/2. (18) 

In this case the Hermitian matrix J~q~ = Jo has the desired property ft = 
ft Jo. This means that the current density is equal to p = ft Jo. In addition, all 
of the gamma matrices now have this property (14) r£ ft = ft Jfe k = x,y, z 
which follows from the equation Jfe = i[bk Jo] (3) and equation (13) ft bu = 
b+ft. 

The characteristic equations and the spectrum are the same for the Hermi- 
tian H-solution and symmetric U-solution, but the corresponding characteristic 
equations for the matrices pj are different and the eigenvalues are not positive 
definite any more [5], that is there are many ghost states. The positive norm 
physical states are lying on the quasilinear trajectories of different slope and 
the negative norm ghost states are also lying on the quasilinear trajectories. 
There are also tachyonic solutions [5], which appear in Majorana equation as 
well [9, 14, 16, 5]. In [5] it was suggested to "protect" equation from tachyonic 
solutions by setting some of the transition amplitudes 7J r to zero. 

In the case when some of the transition amplitudes 7J r in (18) are set to 
zero (we again use "superposition" property of the equation (8) [5]) 

7 1 1 = l] 3 - if 5 - •••• = 7 1 1 - 7 4 3 - 7' 6 = = (19) 

and all other elements of the Jo matrix remain the same as in (18) we have a new 
Zi-solution with the important property that J^ is a diagonal matrix and that 

3 Computing the traces and determinants of these matrices one can get the following 
general relation for the eigenvalues e\ ■ ... ■ <^ +1 / 2 = b e i + ••■ + e j+i/2 ~ j(2j + !)• 
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the antihermitian part of Jfc anticommutes with To. Thus in this case we recover 
the nondiagonal part of the Dirac commutation relations for gamma matrices 
{7o,/\} = k = x,y,z. For the solution (19) one can explicitly compute 
the mass spectrum and the slope of the trajectories [5] 

„,, 2M 2 f - (2n - l)j + n(n - 1) _ „ 

M " = — J -(n -l)/2 " n=1 ' 2 '- (20) 

where j = n + l/2,n + 5/2, The string tension cr„ = l/2ira n varies from one 

trajectory to another and is equal to 

1 2M 2 

27ra„ = — = n=l,2,... (21) 

a„ n 

Thus we have the string equation which has trajectories with different string 
tension and that trajectories with large n are almost "free" because the string 
tension tends to zero. The smallest mass on a given trajectory n has spin j = 
n + 1/2 and decreases as M 2 (j = n + 1/2) = n f^ 3) ■ 

The other solution, ^-solution, which shares the above properties of S\- 
solution is (18) with 

7 12 = 7 34 = =Q 7 i 2 = 7 3 4 = ___ = _ (22) 

The difference between these last two solutions is that in the first case the lower 
spin is j = 3/2 and in the second case it is j = 1/2. The unwanted property 
of all these solutions E, Si and S 2 is that the smallest mass M%(min) tends 
to zero and the spectrum still has an accumulation point. We have to remark 
also that both equations, S\ and E2, which correspond to (19) and to (22) have 
natural constraints [5]. 

Dual equations. The unwanted property of the Z'-solutions, that is the de- 
creasing of the smallest mass on a given trajectory, can be avoided by dual 
transformation of the system [5] (see section "Dual representations"). Under 
the dual transformation = (jo; A) — > (A; j ) = dual the representation = 

(0 I ^,---,0 i ,0 1 ,---,0 N ) is transformed into its dual dual = (A; -5/2) 

(A; —3/2) (A; -1/2) (A; 1/2) (A; 3/2) (A; 5/2).... and we are lead to take A 
to be half-integer and to A = 1/2 in order to have the Dirac representation in- 
corporated in 4 . The solution which is dual to E 2 (18) and (22) is equal to 
[5] 



7 ; +lr = 7r +1 =7r^=7r +1 = V(J-^l) r> j + 3/2 (23) 

where j = 1/2, 3/2, 5/2, r = 2,4, 6, .... and the rest of the elements are equal 
to zero 

7 1 1 = 7 j 2 = 7 | 4 = .... = 7 1 1 - 7 1 4 - 7 S i = = 0. (24) 



These representations do not coinside with the ones in Ramond equation [11]. 
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The Lorentz boost operators b are antihermitian in this case = —bk [5], 
and therefore the J^ matrices are also antihermitian r£ = — J\. The matrix Q 
changes and is now equal to the parity operator P, the relation f2 r+ = Q 
remains valid. The diagonal part of 7^ anticommutes with r as it was before 
{r , r^} = k = x,y, z. The mass spectrum is equal to 

n j + („ + i)/2 v ' 

where n = 1,2,3,.. and enumerates the trajectories. The lowest spin on a given 
trajectory is either 1/2 or 3/2 depending on n: if n is odd then j min = 1/2, 
if n is even then j min = 3/2. This is an essential new property of the dual 
equation because now we have an infinite number of states with a given spin j 
instead of j + 1/2. The string tension is the same as in the dual system (21). 
The lower mass on a given trajectory n is given by the formula (j = 1/2) is 
M 2 (j = 1/2) = 4 f (2 " + ^jf +3) -> (4M) 2 and the spectrum is bounded from 
below by positive mass. 



Generalized wave equation. The last Z'|" o/ -equation has the property that 
only the diagonal matrix elements of the anticommutator {/o,!^} are equal 
to zero < j, m, r\{r , T z }\r, j, m >— 0, and that nondiagonal elements are 
not equal to zero < j - l,m,r\{r , T z }\r, j,m >— isj j 1 — m 2 e r [(7^ 



(7J-1 ) ]■ Let us search the solution of the Majorana commutation relation (8) 
in the same S^^-ia-cohy form (9) but with additional nonvanishing antidiagonal 
matrix elements 7^ r . The solution has the form [6] 



r r _ r r _~r+l r+1 _ _ r+1 r+1 _ 3 + I/ 2 



where j = 1/2, 3/2, 5/2, r = 2,4, 6, .... and r > j + 3/2 and one can check 
directly that To is the solution of (8). These additional matrix elements in Jq will 
not change the diagonal matrix elements of the anticommutator but will cancel 
nondiagonal matrix elements [6] < j — l,m,r\{ro,r z }\r,j,m >= 0. One can 
check this fact also using the relation {Jo r z } = i [b z /g], which follows from 
(3). Using the relations (3) F y — —i \T Z a x ] and [i"b a x] — one can see that 
{In, r y } = and in the same way that {-To, r x } = 0. Finally using the relation 
(3) Jfc = — i [Jo bk] one can prove by direct calculation that {J 1 ^, Jl} = for 
k 7^ I and then using (8) and the fact that [bk -T^] = one can prove that 
r k =~ r o, therefore [6] 

{r^r v } = 2 VliV rl (27) 

where r§ is a diagonal martix. Now the theory is Hermitian in all orders ofv/c. 
The mass spectrum is highly degenerated and is given by the formula 

M] = M 2 (l - J^) r>j + 1/2. (28) 
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New mass terms (a • b) r 5 and (a 2 — b 2 ) can be added into the string equation 
(2) in order to increase the string tension 

{ir^d** M(a-b)r 5 - ff M (a 2 - b 2 ) } <f =0, (29) 

where (a • b) and (a 2 — b 2 ) are the Casimir operators of the Lorentz algebra. 
Now all trajectories acquire a nonzero slope 

4r 2 — 1 

M* = M 2 _^( r -l/2) 2 (l + 2 5 (r-l/2)) 2 r > j + 1/2. (30) 

where r = 2,4,6,...., j = 1/2,3/2,5/2,...., thus M 2 > M 2 (j + l) 2 . Thus the 
spectrum of the theory consists of particles and antiparticles of increasing half- 
integer spin and masses. On every mass level there are particles with spin running 
from j = 1/2 up to j max = r — 1/2, r = 2, 4, .. and the degeneracy is equal to 
d r = 2r — 1 = 2j max . This is in contrast with the exponential degeneracy in case 
of superstrings. The tachyonic solutions which appear in Majorana equation (see 
(20) in [9]) do not show up here. 

We can introduce now the interaction with gauge field using covariant deriva- 
tive 77^ = id^ + gA^ and see that there are no obvious inconsistencies which are 
characteristic to high spin equations in the background field (see W.Pauli and 
M.Fierz in [14]), because the eqaution (2) can be transformed to the form 

{ n 2 r 2 + l -^r^r v - m 2 } <p = o (31) 

and does not produce additional constraints. 

In conclusion I would like to acknowledge Holger Nielsen and Konstantin 
Savvidy for stimulating discussions. I wish to thank Costas Kounnas and George 
Koutsoumbas for support. I would like to acknowledge the organizers of the 
conference for inviting me and for arranging an interesting meeting. 
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